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A  NOTE  ON  INCOMPLETE  INTEGRALS  OF  CYLINDRICAL  FUNCTIONS 


INTRODUCTION 

The  class  of  cylindrical  functions  C  includes  Bessel  functions  of  the  first  kind  J ,  modified 
Bessel  functions  /,  Bessel  functions  of  the  second  kind  or  Neumann  functions  K  (or  A),  Bessel  func¬ 
tions  of  imaginary  argument  or  MacDonald  functions  K ,  and  Bessel  functions  of  the  third  kind  that 
include  Hankel  functions  ot  tne  first  and  second  kind,  H(l)  and  H[1\ 

The  general  incomplete  Lipschitz-Hankel  integral  of  cylindrical  functions  C  v(z)  is  defined  as  the 
function  of  two  complex  variables: 

Cefi(a,z)mfoea‘tllC,(t)dt.  (1) 

Here  the  symbol  e  denotes  the  presence  of  the  exponential  function  and  /*,  v  may  be  complex. 
Analogously,  we  define  integrals  that  contain  the  functions  sin  (at)  and  cos  (at)  in  place  of  exp  (at): 

Cs  a ,  z)  =  Jo  sin  (at)  t*C  „(t)dt  (2) 

Cc  (a ,  z)  =  f  cos  (at)  t*C ..(t)dt .  (3) 

To  assure  convergence  of  Ce  (a  ,  z )  and  Cc  (a  ,  z),  it  is  necessary  that  Re  (n  +  1)  >  |  Re  v  |  when 

C  =  K ,  Y ,  HiX) ,  Ha)\  Re  (1  +  n  +  v)  >  0  when  C  =  /,  J .  When  n  =  p,  we  define,  for  exam¬ 
ple,  C(.  =  Ce  where  for  convergence  Re  n  >  -  1/2  for  all  C. 

Integrals  of  the  type  given  by  Eqs.  (1)  to  (3)  occur  very  often  in  applied  ruailicmatics.  Agrest 
and  Maksimov  [1]  have  found  representations  for  C(,  (a .  z).  Cx  (a,  z),  and  C(.  (a,  z)  using  incomplete 
cylindrical  functions.  In  this  report  we  give  representations  for  C(.  ( a,z ),  Cv  ( a,z ),  and 
C(  (a,  z)  using  only  the  Kampe  de  Feriet  double  hypergeometric  functions  Fiihol-*,  >'!• 

PRELIMINARY  RESULTS  AND  DEFINITIONS 

To  begin,  we  summarize  some  results  that  are  found  in  Ref.  2.  p.  85:  Let  a  and  b  be  arbitrary 
constants, 

F„(c)  -  al.(z)  +  bef"Kr(z) 

G,(c)  =  aJ  v(z)  +  bY„(z) 


Manuscript  approved  February  24, 


/ :  H  =  F 

1:  H  -  G 


1:  H  =  F 
0:  H  =  G 
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Then 


.  -  *>>b 

j  t^H^t^lt  =  e  -  \in  +  V  -  1)cH,,(c).vm  ,  ,,  ,(<vc)  +  (2/3  -  l)«cH(,  |(c ).v#x_  ,.<«c )1 .  (4) 

where  the  Lommcl  functions  s  are  given  hy 

,/•  *  i 


(n  —  v  +  l)(ji  +  p  +  1)  1 


I: 


H  —  v  +  3  /i  +  v  +  3 


(5) 


Now  defining 


f  1:  C  =  /.  AT  (  \:  C  —  K 

*  =  1-1:  C  -  H ,  J ,  Y  11  =  1—1:  C=H,I,J.) 


we  may  deduce  from  Fqs.  (4)  and  (5)  the  result 


i  i 


I yCrU)dt  =  -  \C„(z)  ,F: 


H  —  i>  +  I 


I; 


H  —  i'  +  3  /i  +  v  +  I 


t-- 


+ 


»;cC,,  ,( - ) 


/t  +  i>  +  1 


.F, 


r  +  3  ^  +  i>  +  3 

t  ’  ->  'a 


(6) 


We  define  the  Kampe  de  Feriet  double  hypergeometric  functions  I.  and  (J  and  give  associated 
generating  relations  |3,  4|: 


/- 1«.  it  y,  «;.r.  v|  s  FV0' 


o.o 


_ :  <v;  13: 


a ,  v 


y.  o: 


v  <  oo  .  v  <  oo 


/3.  7 ;  n,  i’,  X;  v,  y  |  =  F^:,1, 


<v,  /3: 


v ,  \' 


>'■ 


A  <  00  ,  |  V 


OO  ( ( X  )  in 

/•In-,  (3;  y,  <S;  v  .  y|  -  V  — -  --  t  y.w  t  (V  v| 


*  (<v)„,((3)m  v"> 

('!<*.  /3,  7;  /i.  »>,  X:  A  ,  y|  =  £  77.7"T..r,~  TTi  |F’I>:"'  +  /<- w  t  e:  V 


0  (/<),„ (id,„(X)„,  hi! 


(7) 


It  is  easy  to  see  that  the  function  /,  is  a  special  case  of  (J: 

(2 1 a,  X,  (3;  7,  X.  X;  v  ,  v  |  =  F  |<r,  (3;  7,  X;  v  ,  y  | . 
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For  brevity  we  define  the  parameu  r  lists 


.  ,  ,  u  +  v  +  t  a  —  v  +  1  ,  a  +  a  +  1  u.  —  i>  +  3  1 

A  ^  „)  =  - - -  ,  - - - ,  1 ;  - - - ,  - - -  ,  — 


.  .  .  fi  +  v  +  1  u  —  v  +  l  .  u  +  v  +  3  a  —  v  +  3  1 

^2(m.  v)  =  - - - , - - - ,  1; - - - , - - - ,  — 


») 


_  fi  +  v  +  2  ji  -  v  +  2  g  +  v  +  2  n  —  v  +  4  3^ 

^  ’  'A  ’  ’  ^ 


b,<„. »)  .  t±±±l .  ,  | .  !L±JL±A  ,  juuuLI  ,  2 

7  2  2  2  ?  2  2 


1  1  3 

^i00  =  —  +  E  y  +  y 


1  3  3 

d2w  =  y  +  m>  i ;  y  +  y 


_  ,  ,  1  +  g  +  a  2  +  g  +  a  2  +  g  +  a  3  +  g  +  a 

£,(>..  ->  -  — 5 — . - 5 -  ■  — 2 — ■  — 2 — • 1  +  " 


F.2(h ,  a)  = 


+  n  +  a  2  +  n  +  a  j  3  +  n  +  a  4  +  n  +  a  |  + 

2  ’  2  ’  ’  2  ’  2  ’  V 


y  +  n,  1;  1  +  n,  y  +  n 


E  2(7*)  =  y  +  l;  2  +  /*,  y  +  m 


REPRESENTATIONS  FOR  C(,  (a,  c),  Cv  (a,  z),  C(  (a.  z) 

Substituting  the  Maclaurin  series  for  exp  («/)  in  Eq.  (1)  and  splitting  into  even  and  odd  terms 
we  obtain  on  integrating  term  by  term 


00  ..1+2  n 


E^rt"ltl-cp(/wr+  E  tt 


„=<>  (1  +2n)! 


Then  using  Eq.  (6)  and  the  generating  relation  Eq.  (7)  we  obtain  after  a  tedious  but  straightforward 
computation  the  principal  result  of  this  note 
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Ce  (a  ,  Z)  =  z'  +  "C Jz) 


H  —  v  +  1 


Q 


.  a~Z  gz 


4  ’  4 


+  — 

/i  —  v  +  2 


*■; 


a2Z2  _g 

4  ’ 


+  T]Z2  ^C,_i(z) 


j  2  2  >  ^ 

_ o  \A  .  “  z  iii 

(/i  +  */ +  l)(/t  -  j/ +  1)  2’  4  ’4 


1 


Since 


+ 


az 


(n  +  i>  +  2){fi  —  v  +  2) 
1 


fl  ■  fl2"2  ^2 


CV .(a’  z)  -  Yi  i  c<v/'«.  z)  -  Ce{-ia,  Z) 


we  may  write 


Cc„.(a'  z)  =  j  j  Ce.  -{ia'z)  +  C*.  z) 


Cs  (a  ,  z )  = 


az 


2  +  i 


H  -  v  +  2 


C„(z)Q 


2  2  y  1 

~a  z  gz 
4  ’  4 


a,; 


rjz 


H  +  i>  +  2 


C„.i(z)Q 


R  - a2z 2  gz2 

4  4 


C(  {a  ,  Z )  = 


-i  +  A* 


^  -  v  +  1 


r 

2  2  v  1 

1  C,iz)Q 

4  a  *•  £' 

/i  1  ,  , 

1  4  4 

+ 


M  +  „  +  l 

=or  n  =  v,  Eqs.  (8)  to  (10)  reduce  to 


C„-,(z)Q 


A,: 


2  2  v  ■> 

-a  z  gz 


4  ’  4 


Q(a.c)  =  c1  4  "C„(z)  U 


2  2  >.  2 

«  Z  gz 


+  —  Q 
2 


41.  ic 

4  4 


+  nz2*^_,(c) 


I  +  2m 


Dz 


2  2  >  *> 
«  Z  gz" 


4  ’  4 


+  — — —  Q 

4(1  +n) 


/*); 


2  2  >  *> 
«Z~  gz- 

4  ’  4 


Cv/«  ’  z)  =  -y  rtz2  4  "  j  C„(z)(? 


i(/l  m); 


2  2  >  ' 
-a  z  gz" 


4  '  4 


+ 


C„  -  ,(z)C 


J  J  y  ’ 

n  (  ,  -a  z  kz 

ft):  — - —  ,  -i— 


in 
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Cc(a.z)-zl^\c^)L  +-M-Cf.l(z)L  D2;  -^1 . 


Defining  J+  =  J ,  J  s  /,  it  is  interesting  to  note  that  we  may  also  write  [6] 


Je±(a,z)  = - —  —  — -  iQ 

2\\  +  g  +  y)r(l  +  v)  t  4  4 


■>  *> 

az  ~  r  T Z “ 

O  L  7,  , 

2  +  g  +  p  ‘  4  4 


z(r/2)V'; 

(i  +  2/x)r(  l  +  n) 


(14) 


Here  the  Bessel  functions  J  *  do  not  appear. 

REDUCTION  FORMULAS  FOR  L  AND  Q 

Many  special  cases  of  Eqs.  (II)  to  (14)  may  be  obtained  in  one  form  or  another,  provided  we 
know  a  reduction  formula  for  either  L  or  Q .  We  summarize  some  known  relevant  reduction  formu¬ 
las  [3-6): 

L[ck,  (3;  7,  6;  c.  =  ,F2[ a  +  j3;  7,  6;  z\ 


2 g  sinh  ;  +  cosh  ; 
1+2  /x  c  1+2  g 


Q  fl2(g.  g);  -y- . 

4  4 


I  +  p  4  L 

TTSTTI  cosh:  - 


-  ro+gj/^c) 


0  b,(m.m);  7.  7 


2  1  I  cosh 

^7  2"  — 


+  sinh  c  -  —  TO  +  g)/H  .  \{z) 


Other  properties  and  reduction  formulas  for  L  and  Q  arc  found  in  Refs.  3-6. 

APPLICATIONS 

Of  interest  in  applications  are  the  functions  Je  (a,z),  /,.  («,;),  Yc  (a.z).  and  AT(.  (zz ,  c ) - 
Je  (a.  z)  and  Ye  (a ,  c)  occur  in  problems  in  the  theory  of  diffraction  in  optical  apparatus  |1.  p.  227). 
The  function  le(a ,  c)  plays  an  important  role  in  the  study  of  oscillating  wings  in  supersonic  flow  and 
arises  in  the  study  ot  resonant  absorption  in  media  with  finite  dimensions  |1,  p.  195).  Kr  (a,  z) 
occurs  when  the  statistical  distribution  of  the  maxima  of  a  random  function  is  applied  to  the  amplitude 
of  a  sine  wave  in  order  to  calculate  the  distribution  of  its  ordinate.  This  latter  distribution  is  of 
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interest  in  the  study  of  the  scattered  coherent  reflected  field  from  the  sea  surface  [7,  8],  Since  the 
functions  Ce  (a ,  z)  are  of  some  importance,  by  using  Eq.  (11)  and  defining 


L{(x,y)  =  L  y,  1;  |;jc,y 


L(fa,y)mL  y,  U~,  |;jc.v 


Q\(x,  y)  =  Q  l,  1,  l;  1,2,  j-.x ,y 


0oU,y)  =  0  1,  1,  1;  2.  2,  y;  .r,  y 


we  obtain 


Kei(u,z)  =  zK0(z)  U 


q2z2  z2 

4  *  4 


7  f  7 

az  a”  7^  7 " 

+  —  0 ,  — —  ,  — 
2  ^  4  4 


//  ^  „  f  z2 

\L°  v  4  ’  4  J  +  4  l  4  '  4 


Veii(a,z)  =  zYq(z)  \l,  +~<2i 


4  4 


+  +  T*«[iTi'¥ 


J'ja.zi-z/jz)  U,  +  f  C.  -J- 


+  c27,(c)  Uo  +f  <?«  £f~>  "7“ 


/,„(«.  ;)  =  zl0(z)  U,  Y  +  yG.  ^T-f 


:2/.fc)  Uo 


7  7  1 

a 

4  '  ~4 


7  ■>  7 

4-  a  z  A_ 

4  4  ’  4 


RS 
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The  equations  for  H}']  and  H™  are  the  same  as  those  for  Ye  or  Je  with  Y  or  J  replaced  by  H{[)  or 
H[1).  Further,  from  Eq.  (14)  we  have 


J  (a,  z)  =  zeaz 


2  ’  2  ’  ’  4  ’  4 


—  L 
2 


1,  — ;  — .  2;  — 
2  2 


l  .  z ) 


L 


Here  we  have  used  the  properties  of  L  that 


L[a ,  |3;  7,  5;  x ,  y]  —  L[  a,  j3;  6,  7;  .r ,  y  |  =  L  |/3,  a;  7,  6;  y ,  .r  J. 


The  latter  results  for  Ct,(a ,  c )  should  prove  useful  in  numerical  computation  of  these  functions. 

SUMMARY 


Representations  for  incomplete  Lipschitz-Hankel  integrals  of  cylindrical  functions  using  only  the 
Kampe'  de  Feriet  functions  in  two  variables  F^kol*'  >1  are  given.  In  addition,  known  relevant 
reduction  formulas  for  these  functions  are  provided. 
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